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ABSTRACT 



Aims. With the availability of galaxy distance indicators in weak lensing surveys, lensing tomography can be harnessed to constrain 
the three-dimensional (3D) matter power spectrum over a range of redshift and physical scale. By combining galaxy-galaxy lensing 
and galaxy clustering, this can be extended to probe the 3D galaxy-matter and galaxy-galaxy power spectrum or, alternatively, galaxy 
biasing parameters. 

, Methods. To achieve this aim, this paper introduces and discusses minimum variance estimators and a more general Bayesian ap- 

. proach to statistically invert a set of noisy tomography two-point correlation functions, measured within a confined opening angle. 

^,1^ ■ Both methods are constructed such that they probe deviations of the power spectrum from a fiducial power spectrum, thereby en- 

abling both a direct comparison of theory and data, and in principle the identification of the physical scale and redshift of deviations. 
By devising a new Monte Carlo technique, we quantify the measurement noise in the correlators for a fiducial survey, and test the 
performance of the inversion techniques. 
^ ■ Results. For a relatively deep 200 deg- survey (z ~ 0.9) with 30 sources per square arcmin, the matter power spectrum can be probed 

in ■ with 3 - 6(T significance on comoving scales 1 < kh^'Mpc < 10 and z ^ 0.3. For 3 lenses per square arcmin, a significant detection 

_ (~ lOcr) of the galaxy-matter power spectrum and galaxy power spectrum is attainable to relatively high redshifts (z ^ 0.8) and over 

d ■ a wider A:-range. Within the Bayesian framework, all three power spectra are easily combined to provide constraints on 3D galaxy 

biasing parameters. Therein, weak priors on the galaxy bias improve constraints on the matter power spectrum. 
Conclusions. A shear tomography analysis of weak-lensing surveys in the near future promises fruitful insights into both the effect 
, of baryons on the nonlinear matter power spectrum at z < 0.3 and galaxy biasing (z ;£ 0.5). However, a proper treatment of the an- 

1^ ■ ticipated systematics, which are not included in the mock analysis but discussed here, is likely to reduce the signal-to-noise ratio in 

■ the analysis such that a robust assessment of the 3D matter power spectrum probably requires a survey area of at least ~ 10^ deg^. To 

^N)^ ■ investigate the matter power spectrum at redshift higher than ~ 0.3, an increase in survey area is mandatory. 

^SJ ' Key words, dark matter - large-scale structure of Universe - gravitational lensing 

o: 

fSl . 1. Introduction in the matter density field pm{x,x) relative to the mean matter 

density Pmix), at a comoving position x, is expressed by the mat- 
In the framework of the standard model of cosmology (e.g., ter power spectrum P^(k,x) (e.g., Martinez & Saar 2002) 
Dodelson 2003), structure in the matter density field evolves , , 

over time under the influence of gravitational collapse from a rel- iSm{k,x)^m{k' - (2^) k')P^{k,x) , (2) 
atively homogeneous state, as for instance observed in the small 

temperature fluctuations of the cosmic microwave background where 

(CMB; Komatsu et al. 2011), to the highly structured matter „ 

density field of today. The main driver of the structure forma- S^(k,x) = ( d^xS^ix^x)^^'^ (3) 

tion is the postulated dark matter, which is the dominant matter J 

component. According to the model, dark matter only interacts ^j-e the amplitudes of the density modes with wave-number \k\ 

by means of either the weak nuclear force or gravity. Moreover, , -.u o /m r(3)/ \ ■ »u ni- j u r j 

/ ^ , . ^ , ^ , ,. , or wave length 2n \k\, o J(.x) is the Dirac delta function, and 

the first galaxies were formed from the subdominant primordial / \ j . uV n v f 

. ° , ., ,. • , • , ••,{■••) denotes the ensemble average over all realisations of the 

baryonic gas component, building stars within the gravitational „ . -.^ijt.u i i ..j -.cu 

. , ° „ , , ' ° . ^^r^^,^^ , matter density field. In the cosmological context, density fields 

potential wells of the dark matter (e.g.. Mo et al. 2010), and » » ii u j • 

evolved later on ^ = ' ^' are statistically homogeneous and isotropic. 

The matter power spectrum contains a wealth of information 
To quantify the state of the structure formation process at a directly related to the physics of the matter components in the 
given cosmic time or a radial comoving distance x, the amplitude Universe as well as the primordial density fluctuations (Peacock 
of the fluctuations 1999, and references therein) and information about possible 

modifications to Einstein's general relativity that may affect the 
growth of the density fluctuations (Uzan & Bernardeau 2001; 
Clifton et al. 201 1). While the matter power spectrum at redshift 
Pjjj(x,x) - Pm(x) ~ ~ 1000 is relatively well-understood theoretically and known 

Smix,x) - p^(x) empiricafly from CMB studies (Komatsu et al. 201 1), the power 
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spectrum in the more recent Universe on quasi- to nonlinear 
scales k > 1 /iMpc ' is less well-known: On the theoretical side, 
the magnitude of the influence of galaxies and baryons seems 
to be unclear but probably significant (Zhan & Knox 2004; Jing 
et al. 2006; Somogyi & Smith 2010; Semboloni et al. 2011; van 
Daalen et al. 2011), although predictions for dark-matter-only 
scenarios are accurate (Peacock & Dodds 1996; Smith et al. 
2003; Heitmann et al. 2010). On the empirical side, detailed 
studies of the power spectrum in this regime rely mostly on trac- 
ers (Tegmark et al. 2004), such as galaxies or absorbers in the 
Lyman-a forest, which have an uncertain relation to the under- 
lying matter density field. This provides a strong motivation to 
directly measure either the spatial matter power spectrum or the 
spatial matter correlation function. 

One purpose of this paper is to propose estimators of the 
matter power spectrum based on the weak gravitational lens- 
ing effect to allow a direct measurement of Pm{k,x), especially 
in the nonlinear regime. Weak gravitational lensing (Schneider 
2006, for a review) is a tool to study the large-scale distribu- 
tion of matter by its effect on the shape of distant galaxy images 
("sources"). Images of galaxies are weakly distorted (sheared) 
by the bending of light passing by intervening matter density 
fluctuations ("lenses"). As the effect is only sensitive to the fluc- 
tuations in the gravitational field of intervening matter, it can 
be employed to probe the matter density field without further 
knowledge of its physical properties and without the usage of 
tracers. This makes it an excellent probe to explore the elusive 
dark matter component or the general composition of matter 
on cosmological scales (see for example Clowe et al. 2006 or 
Tereno et al. 2009). 

As weak lensing is sensitive to all intervening matter, it nat- 
urally measures the projected matter density fluctuations within 
radial cylinders. Known lensing-based estimators of the mat- 
ter power spectrum are therefore essentially two-dimensional 
by either estimating the projected matter fluctuations on the 
sky (Blandford et al. 1991; Miralda-Escude 1991; Kaiser 1992; 
Schneider et al. 1998; Hu & White 2001; Schneider et al. 2002; 
Pen et al. 2002) or translating this into a spatial power spectrum 
at one effective distance (Pen et al. 2003). With the availability 
of distance indicators for the sources, one has started to exploit 
the change in the weak lensing shear yix^) as a function of co- 
moving source distance ;^fs. For a fixed lens distance ;^^(J, this is 
expressed by 

, s A(Xs -Xd) ... 

rCrs) ^ — , (4) 
AOcs) 

with fk(x) being the comoving angular diameter distance. The 
weight of a lens in the total shear signal, including lenses at 
all distances, therefore varies with a varying source distance. In 
a cosmic shear tomography analysis, this is utilised to provide 
improved constraints on cosmological parameters (Hu 2002; 
Massey et al. 2007b; Schrabback et al. 2010; Kitching et al. 
201 1) or to map the three-dimensional (3D) matter density field 
(Hu & Keeton 2002; Massey et al. 2007a; VanderPlas et al. 201 1 ; 
Simon et al. 2012; Leonard et al. 2012) by analysing the lensing 
signal for a series of source distance slices. 

Weak lensing also offers the opportunity to investigate the 
distribution of cosmological objects in relation to the matter den- 
sity field. As in the case of the matter power spectrum, one de- 
fines the clustering power spectrum of galaxies 

(6,{k,x)S,ik',x)) = {2nf6^^\k + k') i^P^{k,x) + (5) 



and the cross-correlation power of the matter and galaxy distri- 
bution 

(6s(k,x)6m(k',x)) = (2n?6'-^\k + k')Ps^{Kx) , (6) 

based on the galaxy number density contrast 6g = ng/fit, - 1, 
where Wg is the number density of galaxies, and h„ their mean 
number density. We note that the definition of Pg{k,x) separates 
the Poisson shot noise contribution, originating from the dis- 
creteness of the galaxies, from the number density fluctuations 
(e.g.. Peacock 1999). For galaxies faithfully tracing the distri- 
bution of matter, all three power spectra Ps{k,x), Psg{k,x)^ ™d 
Pg{k,x) are identical. More generally, however, we have to as- 
sume that galaxies are imperfect tracers, i.e., biased tracers, of 
the matter distribution. Galaxy biasing is commonly expressed 
at the power spectrum level by the biasing functions b{k,x) and 
r{k,x) (e.g., Tegmark & Peebles 1998): 

P6g{k,x) = b{k,xy{k,x)P6{k,x) , (7) 
Pg{k,x) = b\k,x)P6(k,x) ■ (8) 

Measuring galaxy biasing addresses the question of galaxy for- 
mation and evolution as deviations of galaxy clustering from 
(dark) matter clustering - as reflected by the biasing functions 
with b,r 1 - constrain galaxy models (Yoshikawa et al. 2001; 
Weinberg et al. 2004; Springel et al. 2005). 

Non-parametric weak gravitational-lensing methods have 
been proposed (van Waerbeke 1998; Schneider 1998; Pen et al. 
2003) and appHed (Hoekstra et al. 2002; Pen et al. 2003; Simon 
et al. 2007) to assess the galaxy biasing for an effective radial 
distance, but to date without exploiting tomography. The formal- 
ism in this paper is kept general enough to constrain the galaxy 
biasing functions b{k,x) and r{k,x) or, alternatively, the power 
spectra Pgg(k,x) and Pg{k,x) as well. 

To enable the employment of this new tomography tech- 
nique, the appropriate lensing data will soon be available. 
Recent, ongoing, and prospective lensing surveys - such as 
the Canada-France-Hawaii Legacy Survey (CFHTLS'), the 
Panoramic Survey Telescope and Rapid Response System sur- 
veys (Pan-STARRS-), the Dark Energy Survey (DES^), and the 
Kilo-Degree Survey (KIDS'*) - will provide between hundreds 
and thousands of square degrees of shear data endowed with 
galaxy redshifts. With data sets of this magnitude in hand, we 
can expect to extract the 3D matter power spectrum and biasing 
parameters of a selected galaxy population, as discussed in the 
following. 

The paper is structured as follows. Sect. 2 derives a set of 
minimum variance estimators, which, when applied to measure- 
ments of angular correlation functions detailed below, yield con- 
straints on the 3D power spectra Ps{k,x), Psg{k,x), and Pg(k,x)- 
In Sect. 2.6, the estimators are applied to a fiducial survey (Sect. 
2.4) to forecast their performance. The estimators are improved 
within a Bayesian analysis in Sect. 3. Sect. 3.2 applies the 
Bayesian methodology to the fiducial survey. Finally, in Sect. 
4 we discuss the results and draw general conclusions. 

As a fiducial cosmology, this paper uses a ACDM 
model (adiabatic fluctuations) with matter density parame- 
ter Qm = 0.27, of which baryons are given by Qb = 0.046, a 
cosmological constant Q.^ = 1 - flm, and a Hubble parameter 

' http://www.cfht. Hawaii .edu/Science/CFHLS/ 
^ http : //pan- Starrs . if a . hawaii . edu 
^ http : //www . darkenergysurvey . org 
http : //www . astro-wise . org/pro j ects/KIDS/ 
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Hq - h 100 km s"' Mpc"'. Length scales are quoted for h - I, 
while the shape parameter of the dark matter clustering adopts 
h - 0.704. The normalisation of the matter fluctuations within 
a sphere of radius 8/i^'Mpc at redshift zero is assumed to be 
erg =0.81. For the spectral index of the primordial matter power 
spectrum, we use «s = 0.96. 

2. Minimum variance estimators of band power 
spectra 

For extensive recent reviews of weak gravitational lensing and 
the weak lensing formalism, we refer the reader to Van Waerbeke 
& Meflier (2003) or Schneider (2006). 

As we show in the following, fixing the angular diameter 
distance /kOf ), the matter density parameter and the relation 
between comoving radial distance and redshift 



r 



cdz' 
H(z') 



(9) 



within a reference fiducial cosmology (where H(z) is the Hubble 
parameter at redshift z), allows us to find a simple linear rela- 
tion between the 3D power spectra and the angular correlation 
functions as long as lens-lens couplings are negligible. This is 
assumed for the scope of this paper. 

The lens and source catalogues from a galaxy lensing survey 
split into A^source and Aliens subsamplcs, each covering different 
radial distance regimes. Quantities computed from the /th sub- 
sample are referred to in the following by a superindex as in 
"q^'^". Likewise, quantities computed from a pair of subsamples 
i and j are denoted by By tomography, we mean a set 

of angular two-point correlation functions obtained from pairs 
of galaxy subsamples. The pair combinations include galaxies 
from both different and the same redshift bins. 



2. 1 . Matter power spectrum 

Here we start with shear-shear correlations between sources 
(e.g., Schneider 2006, and references therein) 



^^^\m) = (rl"(0)ry'(0')) ± {y'i>{0)j'^'(e')) , 



„0'' 



(10) 



where and y^'^ denote the cross and tangential shear compo- 
nent of sources in the /th subsample, which are both evaluated 
relative to the line \0 = 0' - 6 connecting the two sources in 
directions and d'. For shear fields and positions on the (flat) 
sky, we employ the complex notation, which is widely used in 
the literature. The source catalogue is split into A^source subsam- 
ples with known radial distributions p^\x) that define the prob- 
ability of finding a source within the comoving distance interval 

Using Limber's approximation (Kaiser 1992), the ^+\9) in 
the shear tomography can be related to the 3D matter power 
spectrum as in (e.g., Simon et al. 2004) 



(11) 



4c4 



2n 



X 



where the average lensing efficiency is defined as 

r-Xh 



^X Py (X )- 



A(x') 



(12) 



We have also introduced a(x) as the cosmic scale factor at co- 
moving distance the «th-order Bessel function of the first kind 
Ji,{x), and the comoving radius of the observable Universe ;^'h- 

Moving on to expand the 3D matter power spectrum Pij(k,x) 
as a linear combination of basis functions, we define a 2D grid 
with grid points {k,„,Xn) and m € { 1 , . . . , A^a- + 1 ), n € { 1 , . . . , A^, H- 
1). The grid positions are sorted, < k„,+i and^n <Xii+i- The 
grid points mark the interval limits of a band power spectrum. 
To focus on deviations from a model power spectrum, the band 
powers are defined relative to a fiducial matter power spec- 
trum P^'^{k,x) that describes the expected variation in the power 
within the bands or regions outside the grid 



Ps{k,x) = PT{k,x) 



1 + ^ ^ H„rn{k,x) [fs,nm ' 1] 
/7=1 m-1 



(13) 



where fs^,„„ are factors allowing a deviation from the fiducial 
band powers, and 

H„Uk,x) := { J oth'erifc' '""^ ' ^'''"''""^ ^l^) 

defines a 2D top hat function. Therefore, Pij{k,x) is approxi- 
mated by in total A^aA^^ free parameters fs,mn- The boundaries 
{k[,kNi_+i) and {x\,Xn-+\) confine the range that is allowed to be 
modified in comparison to the fiducial model. According to Eq. 
(13), outside this region one has only Ps{k,x) - Pg'^{k,x)- Inside 
the region, a value of ,„„ = 1 uses either the power expected 
from the fiducial model or a different power amplitude. In this 
paper, we focus in particular on evaluating the coefficients fs^,„„ 
and similar coefficients for other power spectra we introduce be- 
low. We note that by setting P^'^{k,x) = 1 and the grid being suf- 
ficiently large, we approximate the matter power spectrum by 
constant powers within the band with absolute amplitudes fs,,„„. 
Adopting a concrete fiducial model spectrum as a reference may, 
however, be more advantageous as it implements trends in k and 
X within the (broad) bands as expected from theory. If the trends 
(not the absolute amplitudes) are more or less realistic, we can 
expect the factors /a,,,,,, to be only slowly changing functions of 
either k or x- Moreover, the fiducial model makes the recon- 
struction more stable as it specifies the power spectrum within 
a A:-regime that is only poorly constrained owing to the limited 
angular range that is covered by the angular correlation function. 
For the scope of this paper, where fs^mn - I is constant through- 
out, we keep the number of bins with N~ - 5 (linear) and = 10 
(log-bins) relatively small. 

By virtue of this band power approximation, Eq. (11) can be 
cast into the linear form 



where we utilise the basis functions 



87rc402 



(i) (i) k„,+lfk(x)S 

. w'ixw (x) r . , , ^ pfid 

J '^^ :7~; I dssJo,4{s)Pfi 



a(x) 



K,A(x)0 



MxW 



(15) 



(16) 



the Bessel function is used for whereas J4 has to be used 
in the case of X^!_^\ A short description of the numerical evalua- 
tion of this integral and similar integrals to follow can be found 
in Appendix A. The fiducial correlation function ^'Jl^ accounts 
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for any possible invariable contributions from the (A:,;^f)-plane 
that are not covered by the grid 



where 



pM 





{k,x) , if 2 H,„ 

m.n= 1 

, Otherwise 



(17) 



X 



(18) 



unity matrix. For a regular matrix A, the straightforward solu- 
tion is = A ' . For a singular A, on the other hand, a solution 
jc^ - k^b is constrained to be as close as possible to Ajc = b. 

For the problem at hand, the pseudo-inverse can yield a bi- 
ased estimator for the unconstrained or degenerate parts of fs- 
To determine the biased elements in we use as reference 



(23) 



where 1 is a vector with the same number of elements as fs but 
with all elements set to unity; this is the estimator in Eq. (20) 
applied to the data vector ^test = XI + ^fij that is expected from 
the fiducial reference model. The elements in ff that are unity 
are either unbiased or biased. 



For sufficiently large grids as in the mock analysis below, how- 
ever, we find this contribution to be negligible. 

In practise, one obtains a measurement of ^^■''(0/) for a se- 
ries of Z € {1, . . . ,Ng} 0-bins and A^soui-ce(A^source + l)/2 combi- 
nations of source subsamples {if). After arranging the series of 
data points as one compact shear tomography data vector ^ with 
NgNsom-ceiNsoarce + 1) elements, the set of equations in Eq. (15) 
for all index pairs (ij) can be written as 



fid 



(19) 



Here, the band power factors fs.mn, compiled within the vector 
/a, and the matrix elements of X, consisting of the X^l^\6i; m, n), 
are arranged to comply with the structure of ^. The constant 
vector ^fid consists of the fiducial correlation-function values 

For a given noisy ^ from observation, this equation has to be 
inverted with respect to f s- As an estimator for /a, we suggest a 
minimum variance estimator (cf. Zaroubi et al. 1995) 



(20) 



which minimises the residual (||/^ - /a)|p), the average over all 
noise realisations, where || . . . || is the Euclidean norm. Compared 
to the simple inverse X"', the estimator can also cope if the sys- 
tem of linear equations in Eq. (19) is over-determined. In this 
case, redundant information is optimally combined based on the 
noise covariance 



The covariance of the estimator is 



fsfs) - {fs){f\ 



[X'Ni'X]-' 



(21) 



(22) 



2.2. Galaxy-matter cross-correlation power spectrum 

A line of reasoning similar to that in Sect. 2.1 also appUes to 
the galaxy-matter cross-correlation power spectrum Pgg(k,x)- A 
projection of the 3D matter-galaxy power spectrum Pga is mea- 
sured by estimating the mean tangential shear about a sample 
of lenses with a number density contrast 

K,(0) =4^-1 (24) 



on the sky, namely (Bartelmann & Schneider 2001) 



r.(|A0|) + irx(|A0|) 



^g(0')r(«) 



(25) 



where y^iAO) (mean tangential shear) denotes the mean shear 
around a lens at separation AO in a reference frame rotated by 
the angle defined by the polar angle (p of AO :- - 6' , fx is the 
mean cross-shear component (mean cross shear), which cannot 
be generated to lowest order by gravitational lensing and may be 
used as reliable indicator of systematics, and n„{0) denotes the 
lens number density on the sky in the direction of 6 and A^g = 
(«g(0)) the mean lens number density. The mean tangential or 
cross shear do not rely on owing to the statistical isotropy of 
the fields. 

For galaxy samples sliced in redshift with sources from the 

;th sample and lenses from the jth sample, the relation of y['^\ff) 
to the 3D matter-galaxy cross-power spectrum is (e.g., Simon 
et al. 2007) 



3//2q„ 



(26) 



2c2 



2n fk{x)a{x) 



h(x) 



N, Nt 



which can be used to quantify the error and the correlation of er- 
rors in the estimate. Furthermore, the estimator is unbiased when 
is correct. 

The estimator will fail when the matrix product in the square 
brackets of is singular, which happens when the problem is 
ill-conditioned. In this context, ill-conditioned can mean that ei- 
ther fs comprises a redshift or length-scale range that is un- 
constrained by ^, or the impact of different bands on ^ is en- 
tirely degenerate. One usually has to modify the {k„,,Xn) band 
boundaries, the redshift slicing, or the fiducial power spectrum 
P^'^{k,x) accordingly to remove the singularity. To evade the Y^'^\0\ m, n) - 
problem of the undefined estimator, we use a singular-value de- 
composition (SVD; Press et al. 1992) for the inversion of X'N^'X 
in order to find a pseudo-inverse. The pseudo-inverse A^ of a 
matrix A minimises the matrix norm IIAA^ - 1 11, where 1 is the 



As before with the sources, we need to know the radial (proba- 



bility) distribution q)^ (x) for each lens slice. 

If we fix the fiducial cosmological model and adopt for 
Psg{k,x) the previous band power Ansatz in Eq. (13) with free 
scaling parameters fsg.mn relative to a fiducial power spectrum, 
we now obtain 



n-l m=l 

with the basis functions 

3//,?Qn, 



(0) 



^0' 



X 



(27) 



(28) 



dx 



W^'V)gy(y) 
A(x)a(x) 



ds s J2{s) Pf 



KJMB 



Aix)d' 



X 
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and the invariable terms 



with the basis functions 



^'•fid^ ^ " 47rc202 



Jo 



CO 



Mx)a(x) J " \MxW 



This system of equations is cast into a more compact form 

7t = Y/,g+7,fid (30) 

for a series of measurements yt'^Vfl/) and the band power- 
spectrum coefficients fsg,mn appropriately arranged inside the 
vectors 7, (A^aA^iens-^source elements) and /^g (N^N- elements), 
respectively. As in the foregoing section, we place all values 
y^il\{Oi) into a constant offset vector, which is here fi f,^. 

Following from this linear relation, a minimum variance es- 
timator for ftjg is now given by 



where 



(31) 



(32) 



is the inverse covariance of the estimator and Ny is the noise 
covariance of the galaxy-galaxy lensing measurement. We note 
that we have chosen the fiducial power spectrum to be identical 
to the fiducial 3D matter power spectrum in the last section for 
convenience. As lenses are probably biased with respect to the 
matter distribution, we should expect fsg,nm in general to differ 
from unity. By construction, fsg.mn is an average of r(k,x)b{k,x) 
for the lenses over the corresponding band relative to the fiducial 
matter power spectrum. By taking ratios with respect to fg^mn, 
we can, however, estimate r{k,x)b{k,x) with respect to the true 
matter power, as done in the Bayesian analysis below. 



2.3. Galaxy power spectrum 

Finally, we describe our list of minimum variance estimators by 
considering the angular clustering of the lenses (Peebles 1980) 



wW)(|A6»|) = 4"(6l')4 (6*) 



(33) 



which probes the spatial clustering of lenses from the ith relative 
to the yth sample (e.g., Simon et al. 2007) 



In 



Mx) 



X .(34) 



Here, Pg{k,x) represents the lens 3D clustering power spec- 
trum, describing the fluctuations in the lens number density. 
We note that for lens subsamples overlapping in radial dis- 
tance, we would expect a non-zero clustering signal when cross- 
correlating different lens samples, i.e., for ; j. 

Compiling a series of <J'^\0i) measurements within a vector 
<x> of A^eMens(Mens + l)/2 elements, we can again relate the data 
vector to a band power representation of Pg{k,x) 

w = Z/g + fjfid , (35) 

where the elements of oj and Z are determined by 

a^'j\e) ^YjYj ^'''^'^ "^-^s,™" + Wfif (») (36) 



(29) Z^'''\d;m,n) ^ 



1 



27r6i2 



(37) 



. g^^Cy)^'(y) r . , , . pfid 

77r\ dii7o(^)^5 



k,„A(x)e 



A(x)0 



and the invariable terms, to be subtracted from the measured cor- 
relation function 



1 



X 



(38) 



X 



Ax 



flix) 



Jo \A(x)0^l 



Moreover, in analogy to the previous sections, the vector /g now 
compiles the band powers /g,„,„ relative to the fiducial matter 
power spectrum Pg'^(k,x) in the appropriate order, and (o^ is a 
constant off'set vector that is a function of the fiducial cosmolog- 
ical model filled with the values wJ/^'(0/). 

By analogy with Sect. 2. 1, a minimum variance estimator for 
/g is 



/g = [Ng'fz'N-'(t<;-6;fid) , 
where 

n:' = z'Nr.'z 



(39) 



(40) 



and N„ are the inverse estimator covariance and measurement 
noise, respectively. The relative band powers /g,„„, reflect the bias 
parameter b^{k,x) of the lenses averaged over the bandwidth rel- 
ative to the fiducial matter power spectrum. 

2.4. Mock survey 

We now move on to quantifying the accuracy of the estimators in 
Eqs. (20), (31), and (39) for prospective lensing surveys. For this 
purpose, we assume an effective fiducial survey area of 200 deg^ 
in which galaxies, sources, and lenses have a frequency distribu- 
tion p,{z)dz in redshift given by (Baugh & Efstathiou 1994) 



(41) 



where zo = 0.7, a - 2.0, and /? = 1.5, and the mean redshift 
is hence I - 0.9. To account for the impact of redshift er- 
rors, the redshifts in this distribution are assumed to be esti- 
mates, such as photometric redshifts based on the photomet- 
ric filter system of the survey. The estimates are unbiased with 
an assumed Gaussian root mean square (r.m.s.) uncertainty of 
cr, - 0.05(1 + z), which is similar to contemporary accuracy lev- 
els (Hildebrandt et al. 2008). The total mean number density of 
sources is set to n = 30 arcmin"^. The intrinsic shape noise of 
the sources has an rm.s. variance of cr^ = 0.3 for each ellipticity 
component. Hence, the given survey parameters are quite opti- 
mistic in terms of the galaxy number density and survey depth, 
reflecting typical parameters of a presumably space-based sur- 
vey. 

Lenses constitute only ten percent of the total number of 
sources, that are randomly chosen from the total source sam- 
ple. The idea is that in reality we select a subpopulation of 
lens galaxies to be studied, rather than the maximum num- 
ber of galaxies. In the mock data, lenses still, however, have 
the same z-distribution as sources. In this scenario, galaxies 
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Fig. 1. Assumed redshift distribution of galaxies and redshift binning 
for the fiducial lensing survey. Shadowed areas highlight the sample 
selections in the distribution of redshift estimates, and the black lines 
show the true redshift distribution within the samples after including the 
effect of a redshift uncertainty of cr, = 0.05(1 + z). As additional survey 
parameters, we have h = 30 arcmin"^ (sources) and ;i = 3 arcmin"^ 
(lenses), an effective survey area of 200 deg^, and a shape noise variance 
= 0.3. All galaxies are unbiased with respect to the dark matter 
density field. 



are unbiased with respect to the dark matter distribution, i.e., 
b{k,x) - f'{k,x) = 1 for all scales k and radial distances 

For the mock analysis, the galaxy sample is split 
into Mens = A^source = 5 (photomctric) redshift bins 
[0,0.3], [0.3,0.5], [0.5,0.7], [0.7, 1], [1,2] with no overlap. 
The last redshift bin is relatively broad, as we do not expect 
to obtain very accurate redshift estimates for galaxies with 
z e [1,2]. The distribution inside the redshift slices is convolved 
with the Gaussian redshift uncertainty cr. to obtain the true 
redshift distributions. Owing to the uncertainty in the redshift 
estimators, distributions of different samples subsequently 
overlap in the radial direction. The analysis assumes that the 
resulting true distributions are precisely known. The redshift 
binning and distribution of galaxies is illustrated by Fig. 1. 

The grid defining the band boundaries of the three power 
spectra Pii(k,x), Pegik,x), and P„{k,x) is defined by the ranges 
k e [0.01 /iMpc-i,50/iMpc-i] with Nk = 10 log-bins and z e 
[0, 2] for - 5 redshift bins. The redshift bins for the grid are 
identical to the bins of the angular correlation functions. In gen- 
erally, the grid redshift limits could be chosen independently, 
and the number of redshift bins might also differ. For the fidu- 
cial 3D matter power spectrum, we adopt the fitting formula of 
Smith et al. (2003). Matter and galaxies cluster exactly as in the 
reference matter power spectrum, i.e., all /„„ - 1. 

2.5. Measurement noise 

The essential ingredients for quantifying the estimators' accu- 
racy are the noise covariances (shear-shear correlations), Ny 
(galaxy-galaxy lensing), and (lens clustering). This includes 
the source shape noise, lens position shot-noise, and cosmic vari- 
ance of the actual signal. For an extensive discussion of the noise 
terms in ^-t we refer to for example Schneider et al. (2002). 
Towards this goal, we make A^i-eai - 4800 realisations of mock 
shear and lens catalogues in 2 x 2 deg^ fields based on the pre- 
vious survey parameters and the matter/galaxy clustering in our 



fiducial cosmological model; shear fields are Gaussian random 
fields and galaxy number density fields obey log-normal statis- 
tics. The details of the mock catalogue generation are laid out 
in Appendix B. Each mock survey field is processed to estimate 
the cosmic shear correlations ^±^'(0) (Schneider et al. 2002), the 
galaxy-galaxy lensing signal y^l'\6) (e.g., Simon et al. 2007) and 
the lens clustering (J-'-'\0) (Landy & Szalay 1993). The corre- 
lators are hence estimated exactly as in a real analysis; refer- 
ences for the employed estimators are given inside the previous 
brackets. The correlation functions are estimated between angu- 
lar separations of € [3", 2 deg] using Ng - 30 log-bins. The to- 
mography measurements for all three correlators are separately 
combined into a data vector rf, for each realisation / = 1 . . . Meai- 
The data vectors have the following structure: 

- ^±''(0): The first half of di contains estimates of whereas 
the second half contains the f_ measurements. Both the ^+ 
and f_ values are sorted as tomography blocks of constant 
index {ij) and ascending index /, / < j. Within a block, values 
are arranged in order of increasing angular separation bin 9i. 
Since there are Ng angular bins, A^source(A^souice + l)/2 tomog- 
raphy blocks, and two correlation functions, the data vector 
realisations di have in total A^flA^source(A^soui-ce + 1) = 900 ele- 
ments. 

- y^i^\0)'- Likewise, GGL measurements are sorted in order of 
tomography block {ij) index with ascending lens sample in- 
dex / and source sample index j. All {ij) combinations are 
allowed, i.e., MensA^source combinations. Inside a block, val- 
ues are ordered by increasing lens-source angular separation. 
This yields in total A^eMensA^soui-ce = 750 elements. 

- (J'^\9): Measurements of the angular lens clustering are 
again sorted in order of tomography bin index {ij) with 
; < j and /, j being the lens sample indices; there are in total 
Mens (Mens + l)/2 block Combinations. Inside a block, mea- 
surements are sorted in order of increasing lens-lens separa- 
tion. In total, this yields MeMiens(Mens + l)/2 = 450 elements. 

Noise-free data vectors with fiducial survey parameters are de- 
picted in Fig. 2. 

As an estimator of the noise covariance based on the mock 
data realisations, we devise the field-to-field variance 

Cest = V {di - d) {di - dj , (42) 

^real t -^^ 

where 

d^ — Vdi (43) 

M-eal ^ 

is the ensemble mean. The resulting matrix Cest estimates the 
noise covariance N^, Ny, or N^^ in a A4 = 2x2 deg^ survey de- 
pending on whether di consists of ^l^\ff), -y|''''(0), or (J''^{6), re- 
spectively. We scale this to the expected noise in a larger A200 = 
200 deg^ survey by multiplying Cest with A4/A200 = 1/50. This 
mimics a survey of area A200 consisting of 50 statistically inde- 
pendent patches of size A4 each. Finally, since the band power 
estimators require the inverse noise covariance, we utilise the 
inverse covariance estimator of Hartlap et al. (2007), which is 
based upon Ces'- 

For the covariance estimation and the mock analysis, lenses 
are a randomly chosen subsample of the sources (ten percent). 
In particular, sources are clustered as the projected dark matter 
distribution up to the level of the Poisson shot noise. The choice 
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Fig. 2. Noise-free data vectors ^ (left), yt (middle) and oj (right). Used are A'sourcc = 5 source redshift slices and N\^„i = 5 lens redshift slices; i,j 
are indices of slices sorted by ascending mean redshift. The data vector elements are sorted by tomography index (ij) and then angular separation, 
Ng = 30. A spike corresponds to one (if) tomography block with varying 6. 



of identical redshift bins for lenses and sources is due purely 
to the means by which the mock survey galaxy catalogues are 
generated. In the estimator formalism, lenses and sources may 
have different radial distributions and the number of subsamples 
may differ as well. The z-slices in the mock data overlap due 
to the adopted redshift uncertainty. Therefore, we have a non- 
vanishing for / ^ j, which can be seen in the right panel 
of Fig. 2. 

2.6. Forecasts 

Using the fiducial survey and the predicted measurement noise, 
constraints on the three power spectra are shown in Fig. 3. The 
figures are obtained by applying the estimators to a noise-free 
data vector (Fig. 2), whereas the error bars are based on esti- 
mates of the measurement noise in the mock survey, that is ei- 
ther Eq. (22), (32), or (40). The data points or the centre of the 
shadowed region therefore simulate the statistical average of a 
band power estimation process for an infinite number of simi- 
lar surveys. The error bars or regions, however, are the predicted 
Icr uncertainties in a typical single 200 deg^ survey. As can be 
seen, the average data points lie exactly on = 1, indicating 
that they were evaluated with an unbiased estimator. 

In general, the constraints involving cosmic shear become 
worse with increasing redshift owing to the decreasing number 
of background sources at higher redshifts. For the matter power 
fs,mn, we find only upper limits in the redshift bins z > 0.5 or 
for scales k < 0.3 /iMpc"', but reasonable constraints for z < 0.4 
andO.5 < fc/i 'Mpc < 10. A similar conclusion can be drawn for 
fss.mn, albeit the constraints are here considerably improved, es- 
pecially at higher redshift. Under our assumptions, the galaxy 
band power spectrum of the lenses, /g,m„, is extremely well- 
recovered. The detection fidelity depends mainly on the num- 
ber of lenses within a redshift bin. In reality, this number is 
likely to be even smaller as in the mock survey, depending on 
the galaxy population selected for investigation (see, e.g., Zehavi 
et al. 2011). 

Since the correlation functions cover only a certain angular 
range, we may anticipate obtaining biased estimates for large- 
scale modes k. This bias can be quantified by utilising both Eq. 
(23) for f'^^ and equivalent bias indicators for the two other band 

powers fg„ and /g. For the band limits and fiducial power spec- 
trum chosen, the bias is smaller than 0. 1 percent and therefore 
no concern to us here. We note, however, that a bias clearly be- 



comes visible for Pf(k,x) = 1 and yt < 0.2/jMpc-'. Therefore, 
a realistic matter power spectrum when used as a reference aids 
the stability of the minimum variance estimators. 

The correlation of the estimator uncertainties are depicted in 
Fig. 4. Uncertainties in the estimates of the relative band pow- 
ers fgjnn clearly become increasingly correlated at higher red- 
shift. Most independent estimates are located in the lowest red- 
shift bin, where the constraints are also tightest (Fig. 3). For 
z > 0.4, the errors evidently become strongly correlated because 
all shear-shear correlation functions are sensitive to the same 
matter density fluctuations in front of the sources: For constant 
k, fs.mn adjacent in redshift are anti-correlated. For neighbouring 
bands in the same redshift bin, errors are also strongly correlated 
for z > 0.4. Therefore, there is little statistically independent in- 
formation in the fg^,,,,, at higher redshift. The correlation of errors 
becomes small for the galaxy-matter band power fsg.nm and is 
almost absent for the galaxy band power /g,,„„. However, all es- 
timates still exhibit strong correlations between adjacent A:-bins 
in the large-scale regime (small k) because the bins are associ- 
ated with density modes that are barely touched by the angular 
correlation functions of limited angular range. 



3. Bayesian analysis 

The minimum variance estimators of the foregoing section are 
quick and easy to apply to the data. Moreover, their statisti- 
cal errors are equally straightforward to assess. However, they 
have the disadvantage that band power estimates can also be- 
come negative for the auto-correlation power spectra Pg(k,x) 
and Pg(k,x), which by their definitions is not allowed. One fre- 
quently finds band power estimates oscillating about zero where 
errors are large. Therefore, the question arises of whether one 
can refine the analysis to explicitly incorporate non-negative 
powers. This refinement may also help to improve the recovery 
of the matter power spectrum, which in the foregoing minimum 
variance analysis was only successful at low redshifts, where the 
majority of band power factors ,„„ had only upper limits. We 
note that the cross-correlation power Psg{k,x) can in principle 
become negative, if the galaxy correlation parameter r{k,x) is 
negative. 
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Fig. 3. A simulated reconstruction of the 3D band matter power spectrum Psik,x) (left), galaxy-matter cross-correlation power spectrum PsB(k,x) 
(middle), and galaxy power spectrum P„(k,x) (right) utilising the minimum variance estimators in Sect. 2. Shown are the average estimated band 
powers /„,„ relative to the fiducial band power spectrum. The shadowed area in the background indicates the 68% confidence region of the lowest 
redshift bin z = 0.15. The data points denote estimates at higher redshifts. The en'or regions indicate the Icr uncertainties of the estimators via Eq. 
(22), (32), or (40). Data points in blue colour or with asymmetric error bars indicate estimates that are consistent with zero. Note that data points 
within the same k-hin are slightly shifted with respect to each other along the fc-axis for clarity. 
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Fig. 4. Correlation matrix of uncertainties in the relative band powers (minimum variance estimators). Roman numbers indicate the red- 
shift bins (z = 0.15,0.4,0.6,0.85,1.5 for I-V, respectively). Within blocks, highlighted by dashed lines, the correlation between A:-bins 
of either identical z-bins (diagonal blocks) or different z-bins (off-diagonal blocks) is plotted. The mean wave numbers of the i-bins are 
k = 0.01, 0.02, 0.05, 0. 13, 0.30, 0.71, 1.66, 3.88, 9. 10, 21.33 in units of h Mpc"' . 



3.1. Method 

A refinement can be achieved in the framework of a Bayesian 
analysis (cf. MacKay 2003), where we determine the posterior 
likelihood 

p(m\dj cc £{d, fw) X p^riox[m) (44) 

of our model parameters m - for instance, the band power values 
fs in the case of Ps{k,x) - for the observed tomography correla- 
tion functions d, which is ^ for the example mentioned; X.{d, m) 
is the data likelihood function and Ppiioi(ffi) the a-priori informa- 
tion on the band power (prior). 



In the case of fs, we solely assume a prior that enforces pos- 
itive or zero band powers 

Ppno.{fs)^Y]Y]H{fs,,„„), (45) 

n=l H7=l 

where H(x) is the Heaviside function of x. A conceivable modifi- 
cation of the prior could consist of including available empirical 
information on the matter power spectrum. For the UkeUhood 
function, we presume Gaussian noise 

(^, /,) (X exp |-i - Xfs - ^m]' N^-' - X/,- - f fid]) ,(46) 
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with a constant noise covariance that does not depend on the 
matter power spectrum. The likelihood functions £.y and i]„ for 
/ag and /g, respectively, are defined accordingly. Furthermore, 
for /g „,„ we employ a prior similar to the /j.,,,,, prior, but allow 
/^E./wn to be negative (negative correlations). 

Another convenience of the Bayesian approach is that one 
can easily accommodate a transformation of model parameters 
in the analysis, yielding the statistics of an alternative parameter 
set. We take advantage of this by also determining the posterior 
likelihoods of band galaxy biasing parameters r,„„) by ex- 
pressing /^g and /g as 



/g,mn — ^m„ fsjnn i 



(47) 



for all m - I . . .Nk and « = 1 . . . A^;. Here the prior ensures that 
the bias factors b,,,,, and the relative band powers fs^,„ are pos- 
itive, thereby directly constraining b{k,x) and r(k,x), provided 
that we employ the same reference fiducial power spectrum for 
the matter, galaxy-matter, and galaxy power spectrum. In this 
variant, we combine the information on all {fs.mn, fsg.nm, fg.mn) to 
constrain the biasing parameters. Consequently, we now have a 
combined posterior likelihood instead of the three separate ones 
in the previous case 



p(^fs, b, r\dj oc ;5pi-ior(r, b,fsj X 



(48) 



We assume that noise correlations between y,, and u estima- 
tors are negligible, hence likelihood function factors are com- 
bined as shown, where {b, r) denotes the biasing parameters and 
(brfs, b^f 5) provide an abbreviated version of Eq. (47). 

We sample the posterior likelihoods numerically by devis- 
ing the Markov chain Monte Carlo (MCMC) technique (e.g., 
MacKay 2003), where the noise covariances of the minimum 
variance estimators / are used as yardsticks to set up the pro- 
posal function for the Metropolis-Hastings algorithm. 

3.2. Forecasts 

Fig. 5 shows the constraints on the fiducial survey (Sect 2.4) 
obtained via the Bayesian analysis. We show only the variant 
that uses biasing parameters to express the galaxy-matter and 
galaxy power spectrum. As before with the minimum variance 
estimators, the constraints are obtained by using a noise-free 
data vector as input data, as shown in Fig. 2. The likelihood 
functions Xf.y.oj assume a level of noise that is similar to that 
m the 200 deg^ fiducial survey. This combination basically sim- 
ulates the typical posterior of a 200 deg^ survey. We plot the 
means of the marginalised one-dimensional (ID) posteriors of 
{fs,nm^b„,„,r,„„) and the rm.s. variance in their MCMC values 
above (upper error bars) or below the mean (lower error bars). 
This defines a 68 percent credibility region for every parameter. 
The prior is defined to constrain the galaxy biasing parameters a 
priori for b{k,x) 6 [0, 4] and r{k,x) 6 [-2, 2] (flat priors). These 
priors broadly delineate the regime that is expected for galaxy 
bias (e.g. Guzik & Seljak 2001; Seljak et al. 2005). 

We note that owing to the subtraction of Poisson shot noise in 
the definition of Pf,{k,x) (Eq. 5), the modulus of the correlation 
coefficient r can be larger than unity. A value of \r\ > 1 indicates 
a discrete galaxy sampling that differs from a Poisson process 
as, for instance, predicted on small scales in the halo model (cf. 
Guzik & Seljak 2001). The Monte Carlo mock data employed 



in this study uses an explicit Poisson process to generate galaxy 
positions out of a smooth number density field (Appendix). 

For fs, a flat prior on a logarithmic-scale was used. However, 
we found that a prior flat on a linear scale yields similar results 
for in our case. This result, which is based on a noise-free data 
vector, shows that the mean of the ID posterior probability dis- 
tribution function (p.d.f.) is not an unbiased estimator of any /,„„, 
but the mean is biased to too low values for ,„„ and biased too 
high for b„,„ and for most factors r,„„. This result also depends 
on the adopted priors, particularly when the constraints from the 
data are weak. For example, the flat prior for b,„„ tends to pre- 
dict values of ^,„„ of around two. We note, however, that here 
the true f„,„ = 1 always falls into the above defined 68 percent 
credibility region such that the bias appears to be at least smaller 
than the rm.s. uncertainty caused by measurement noise. This 
can be explained in the following way. First, the posterior mean 
{film) is an optimal estimator f^j^l in the sense that it minimises 
the average estimator error 



2 ^est 



:= j df,„„ p{f,„„\d) {f::^-f,nnf 



as to the true /,„„ with posterior p{fi„n\d) because 



■est 
mn 



for /;,^' = (/;„„) = d/,„„ p{f,„„\d) fr, 



(49) 



(50) 



Second, the average error when adopting = ifmn) is simply 
the rm.s. variance in the posterior, as can be seen by substituting 
</™>for/;r,:,'inEq. (49). 

As with the minimum variance estimators, the strongest 
constraints are found for the regime 0.5 < fc/i"' Mpc < 10. The 
constraints on the band matter power spectrum, for which we 
had upper limits only for the minimum variance estimates, are 
clearly improved by the assertion of positive auto-power spectra. 
We note the apparent strong improvement for k < 0.5 /zMpc"', 
which is merely a combined effect of the priors on /^_„,„ and the 
bias factor and, possibly, also the information on the galaxy clus- 
tering on large scales from the data. 

In addition, the MCMC results are devised to estimate the 
correlation of the errors as a correlation matrix in Fig. 6. This 
also shows an improvement for fs^,„„. The correlation of errors is 
softened, hence much less pronounced than in Fig. 4, although 
we can still make out the anti-correlation among the bins that 
are adjacent in redshift. The (band) galaxy bias and correlation 
factors b„,„ and r„,„ have errors that are very similar to those of 
fsjnn since the correlation pattern of the latter is imprinted in the 
galaxy biasing parameters owing to the mixing of fs^,„n, fsg,mn, 
and /g„,„ within the Bayesian analysis. 

4. Discussion and conclusions 

This paper has presented methods for constraining from tomog- 
raphy data the deprojected power spectra of the matter cluster- 
ing and the related clustering of galaxies. Moreover, a fiducial 
survey has been employed to determine the prospects for the ap- 
plication of this method to real data. 

The prerequisite for a successful application of this tech- 
nique is the availability of radial (redshift) information for sur- 
vey galaxies, which will be provided by ongoing lensing sur- 
veys and surveys in the near future. On the basis of this radial 
information, sources and lenses can be divided into subsam- 
ples with known radial distributions, which are disjunct slices in 
the optimal case, thereby probing matter and galaxy clustering 
at different effective redshifts and on different physical scales. 
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Fig. 5. Bayesian MCMC analysis of the fiducial survey estimating the matter band power (left; relative to the fiducial power spectrum) and the 
galaxy biasing parameters (middle: bias factor; right: correlation factor). The shadowed area in the background highlights a 68% credibility region 
about the marginalised posterior mean for the lowest redshift bin 5 = 0.15, whereas data points delineate the 68% credibility regions about the 
mean at higher redshifts. The highest, most poorly constrained redshift bin is not included here. The results are subject to some numerical noise 
caused by the adopted MCMC technique. 
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Fig. 6. Correlation matrix of band power and galaxy biasing parameter estimates in the Bayesian framework (see Fig. 5). See caption of Fig. 4 for 
more information. 



From the cosmic shear tomography alone, the spatial matter 
power spectrum can be measured, which represents an advance 
on older lensing methodologies (Schneider et al. 2002; Pen et al. 
2003; Bacon et al. 2005). Augmenting this with measurements 
of galaxy-galaxy lensing and angular lens clustering for all pair 
combinations of galaxy subsamples, constraints on the galaxy 
biasing can also be made as a function of redshift and length- 
scale. This improves older non-parametric lensing techniques for 
probing galaxy biasing (Schneider 1998; van Waerbeke 1998) 
since tomography correlations can be directly traced back to the 
underlying 3D power spectrum. In particular, the calibrations of 
bias parameters due to projection effects developed by Hoekstra 
et al. (2002) are no longer required. This promises to provide a 
simplified comparison to theoretical galaxy or dark matter mod- 



els (cf. Yoshikawa et al. 2001; Weinberg et al. 2004; Springel 
et al. 2005). 



4.1. Methodology 

We have shown in this paper that by approximating the 3D power 
spectra by band power spectra relative to a fiducial band power 
spectrum, there is a straightforward linear relation between the 
observable angular correlation functions (shear-shear correla- 
tions, mean tangential shear about lenses, and lens clustering) 
and the 3D power spectra (matter power, matter-galaxy cross 
power, and galaxy power), where Limber's approximation ap- 
plies. This study proposed two approaches to inverting this Unear 
relation in a statistical way: 
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- The first approach makes use of minimum variance estima- 
tors that provide estimates of the relative 3D band powers. 
The approach is similar to that of Pen et al. (2003), but now 
permits the full amount of tomography information to be 
used. By doing so we can now, therefore, observe the time- 
evolving power spectra rather than the average power at an 
effective depth of the survey. Furthermore, by choosing a the- 
oretical power spectrum as reference, the minimum variance 
estimators allow us to compare a theoretical band power di- 
rectly to the data. Differences between the theory and data 
are thereby revealed by a relative band power significantly 

- The second approach is a full Bayesian analysis for which 
we can, in addition, express the galaxy-matter and galaxy 
power spectrum in terms of the linear galaxy biasing func- 
tions h{k,x) and r{k,x)- The biasing functions are related to 
the angular correlation functions in a nonlinear fashion. The 
Bayesian framework automatically accounts for this nonlin- 
earity in the posterior likelihood of the biasing parameters, 
being the output of the analysis. Moreover, we can easily 
implement the prior knowledge that the matter and galaxy 
power, by definition variances, have to be positive or zero. 
On the other hand, the disadvantage of the Bayesian ap- 
proach is that we have to specify the complete noise statis- 
tics of the measurements, whereas the minimum variance 
estimators formally only require the noise covariances. For 
the scope of this paper, we adopted a multivariate Gaussian 
noise model. This does not exactly hold true (Schneider & 
Hartlap 2009; Keitel & Schneider 201 1), but a more sophis- 
ticated noise model can be built into the analysis by modi- 
fying the likelihood function in Eq. (46) in future applica- 
tions. Moreover, the Bayesian framework also allows us to 
add prior information about the matter power or galaxy bias- 
ing from other observations. 

Both approaches rely on a known angular diameter distance, 
f\Sx)^ in addition to a specified relation between radial comov- 
ing distance;\^ and redshift for a given fiducial cosmology model. 
Although this may sound like a strong restriction, those geomet- 
ric properties of the Universe are known to great accuracy and 
are consistently inferred by a wide range of observations (e.g., 
Komatsu et al. 20 11). The nonlinear structure of matter at z < 1, 
on the other hand, is prone to a larger uncertainty as discussed, 
for instance, by van Daalen et al. (201 1). 

4.2. Mock survey 

To estimate the performance of the estimators, we assumed a 
fiducial survey with a 200 deg^ area, a mean galaxy redshift of 
z - 0.9, reaching as deep as z = 2.0 for usable redshift esti- 
mates, and a source number density of h = 30 arcmin"- (Sect. 
2.4 for more details). Lenses are a subsample of ten percent of 
the source catalogue, but have an identical z-distribution. This is 
because a smaller galaxy population is usually studied as lenses. 
We assumed that the angular correlation functions can safely 
be measured between 3 arcsec and 2 degree. Separations much 
smaller than 3 arcsec will be difficult to measure owing to the 
apparent size of smeared galaxy images or the pixel size in the 
instruments. Separations larger than ~ 2 degree are principally 
possible, but this would extend into a regime where the flat-sky 
approximation, used for the formalism of this paper, would begin 
to fail. This may be remedied by a full-sky treatment comparable 
to Castro et al. (2005). 



The fiducial parameters were chosen to be somewhere in- 
between typical figures of current or near future surveys. The 
assumed survey depth is, however, probably more on the side of 
deeper space-based lensing surveys such as COSMOS^ or the 
upcoming Euclid mission*", rather than on the side of shallower 
ground-based surveys such as CFHTLS, KiDS, or Pan-STARRS. 
On the other hand, the adopted number density of sources is 
smaller than expected for a spaced-based mission but still higher 
than the typical number density achieved in ground-based sur- 
veys (~ 10 arcmin"^). Moreover, the modest fiducial survey area 
of 200 deg^ reflects the order-of-magnitude for a contemporary 
ground-based lensing survey, that is capable of delivering cosmic 
shear data such as the CFHTLS. In the not too distant future, this 
figure is likely to be increased by another order of magnitude. 

4.3. Simulated noise covariance 

The predicted noise covariance of the correlation functions was 
derived from realisations of mock surveys (Appendix B). In 
analogy to the Monte Carlo method put forward in Simon et al. 
(2004), I included a new feature to account for the clustering of 
galaxies, obeying a log-normal statistics, and to include a mock 
galaxy-galaxy lensing signal. This approach is imperfect in sev- 
eral respects; 

- The assumption of Gaussian statistics for the cosmic shear 
fields is inaccurate for nonlinear scales (roughly ten ar- 
cmin or smaller), which leads to an underestimation of the 
shear noise covariance (Semboloni et al. 2007; Hartlap et al. 
2009; Hilbert et al. 201 1). The adopted log-normal statistics 
(Bouchet et al. 1993) of the galaxy clustering, on the other 
hand, is presumably a good approximation of the galaxy- 
galaxy lensing and galaxy clustering covariance. 

- As realisations on grids are performed (grid sizes are 4 x 
4 deg^ from which the smaller patches are cropped), we miss 
angular modes on scales larger than > 4 deg (Simon et al. 
2004). Computation time constraints in our approach do not 
allow much larger fields with more galaxies. 

- The covariance is estimated on a patch-by-patch basis, which 
assumes a discontinuous, patchy survey rather than a con- 
tiguous survey area where galaxy pairs across patches could 
also be used to estimate the angular correlation functions. 
This reduces the effective number of galaxy pairs and yields 
more noise than for a realistic contiguous survey area. 

Although the last effect presumably partly compensates for the 
two others, we should expect an overly optimistic noise covari- 
ance and hence possibly overly optimistic constraints on the 
band power spectra. In the future, a more accurate estimate 
could be obtained by ray-tracing through cosmological structure 
formation simulations (White & Hu 2000; Hilbert et al. 2009; 
Kiessling et al. 2011; Li et al. 2011). As the noise covariances 
for the tomography correlators are necessarily large in size (Fig. 
2 for bin numbers), however, a large number of statistically in- 
dependent realisations would be required (this paper uses 4800 
realisations, thus in total 19200 deg^) to properly estimate the 
inverse noise covariance (Hartlap et al. 2007). 



http : // cosmos . astro . caltech . edu/ astronomer/hst . 
html 

* http://sci.esa.int/euclid, see also Laureijs et al. (2011) 
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Fig. 7. Predicted signal-to-noise ratio of the relative band powers f„,„ 
for different redshift bins (see key for line styles in the top panel). 
Minimum variance estimators and Bayesian analysis give comparable 
results except for the matter band power /j„,„ (top panel), which here 
combines constraints from lensing and galaxy clustering with priors 
b <A and |r| < 2 in the Bayesian approach. The bottom panels shows 
the signal-to-noise of the minimum variance estimates. 
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Fig. 8. Signal-to-noise ratio per band of the reconstructed matter power 
spectrum /a„„, in a Bayesian analysis, either for an analysis of the shear 
tomography alone ("independent"), or when combining galaxy cluster- 
ing, lensing with priors < 4, |r| < 2 ("combined"). Plotted are the three 
first redshift bins z = 0. 15, 0.4, 0.6 (top to bottom line). 



4.4. Reconstructed band powers 

The fiducial survey predictions are summarised by Figs. 3 and 5 
for the minimum variance and Bayesian approach, respectively. 
The latter does not include the highest redshift bin, z - 1.5, be- 
cause of its poor constraints. Figs. 4 and 6 visualise the cor- 
relation of errors in the estimates. The correlation of errors 
is quite strong for /^,,„„, but considerably weaker for fsg,mn 
and /g „,„ (similar between the minimum variance and Bayesian 
technique). Compared to the minimum variance estimator, the 
Bayesian analysis eases the strong coiTelation of errors in the 
matter band power spectrum,. 

We adopted A^, - 5 redshift bins and A^^^ = 10 k-bins on a 
logarithmic-scale. The MCMC algorithm is well-suited to fitting 



many free model parameters simultaneously (here A^a-A^; = 50), 
but converges slower if their number becomes larger. For this 
reason, a small number of bins is desirable. Because adjacent 
bands are already moderately to strongly correlated in the k- or 
z-direction, especially for fsj„„, a much finer binning than the 
one adopted is probably not necessary, although this should be 
quantified in an objective way for future applications. A princi- 
pal component analysis (e.g., Tegmark et al. 1997), based on the 
minimum variance estimator noise, could be employed to iden- 
tify the most significant modes in the (fc, z)-plane and find an 
optimised representation of the power spectra for the Bayesian 
analysis. Two interesting binning extremes are conceivable: 

- On the one hand, we could use Nk = 1 and a larger number 
of redshift bins. Using a fiducial 3D power spectrum with 
a reasonably realistic shape but normalised at all redshifts 
to the variance in the local Universe, we could focus on the 
overall growth of the fluctuations from the tomography data 
alone. This would result in an analysis similar to the one in 
Bacon et al. (2005). 

- On the other hand, we could use = 1 and a larger num- 
ber of A:-bins. If the time evolution of the structure growth is 
realistically built into the fiducial power spectrum, then the 
methodology from this paper would infer the deviation from 
the fiducial power as a function of scale averaged over the 
full redshift range. This would focus on deviations from the 
shape of the theoretical power spectrum, taking the expected 
structure growth out of the equation. 

As long as estimates are very noisy - auto-coiTelation band 
powers are consistent with zero in the minimum variance ap- 
proach - the Bayesian approach provides tighter constraints 
(Fig. 7). Both approaches are otherwise comparable. The signal- 
to-noise ratio in the Bayesian analysis equals the mean of the 
marginalised posterior divided by the variance about the mean. 
On large scales A: < 1 /jMpc"', however, the signal-to-noise does 
not fall off as quickly as for the minimum variance estimators: 
it forms instead a plateau. By performing a Bayesian analysis 
for the matter power spectrum combined with galaxy cluster- 
ing and one analysis based on shear tomography alone, we have 
found that the plateau is the result of the asserted positive power 
Peik,x) ^ ("independent" lines in Fig. 8) rather than the effect 
of adding galaxy clustering information and the upper limits to 
b and r ("combined"). 

We expect, on average, the signal-to-noise to decrease with 
redshift, although there are more complicated trends on large 
scales (small k) for /g „,„ because galaxy samples at higher red- 
shift probe larger scales and their galaxy numbers increase first 
with redshift followed by a decline beyond z ~ 0.9. We have 
found that the most significantly measured quantities are the 
galaxy clustering, then the galaxy-matter cross-correlation, and 
then the matter clustering. In reality, we should expect the sig- 
nificance of the first two to go down, however, if an even smaller 
population of lenses is studied than in the fiducial survey. 

The signal-to-noise of the relative band power /j_„,„ drops 
below ~ 2cr on all scales beyond z > 0.4 so that we should 
not expect strong constraints beyond that redshift. For a sur- 
vey shallower than the fiducial survey, this limit is bound to 
decrease yet further The tightest constraints are obtained for 
z ~ 0.2 and about k ~ 2/1 Mpc"' (~ 6cr detection), away from 
which the detection fidelity drops relatively quickly below 3cr 
for individual bands. In conclusion, for z < 0.8 and comoving 
0.5 < kli^^ Mpc < 10 the 3D band power spectra /^g_,„„ and/g_„,„ 
can be recovered most effectively, whereas a significant detec- 
tion of fsjnn is restricted to lower redshifts z < 0.3 inside a range 
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Fig. 9. Predicted signal-to-noise of the galaxy biasing parameters 
(bi,m,r„m) for different redshift bins z = 0.15,0.4,0.6,0.85 (from top 
to bottom line). 



1 < kh-^Mpc < 10. This is in the heart of the regime discussed 
in van Daalen et al. (201 1), hence an application of lensing to- 
mography promises tom improve our knowledge of the effect of 
baryons on the matter power spectrum. 

As we expect the noise covariances to be overly optimistic, 
a robust study of the matter power spectrum would, however, 
probably require coverage of a larger survey area than in our 
fiducial survey. We can hope for a significant signal enhance- 
ment, if the survey area is increased to ~ 10^* deg- or more, as 
expected from full-sky surveys. The signal-to-noise should then 
roughly be boosted by VTO ~ 3.3 for all band power estimates. 
Moreover, boosting the number density of sources, for example 
through space-based surveys, also increases the signal-to-noise 
as shot-noise on small scales is oc H"'^^ (Schneider et al. 2002). 

4.5. Reconstructed galaxy biasing parameters 

Combining the band power estimates, the Bayesian analysis pro- 
vides constraints on the galaxy biasing functions (see the two 
right-hand panels in Fig. 5). Here, the highest accuracy (higher 
than ~ 3cr) is achieved forz < 0.5 and 1 < kli^^ Mpc < 10 (Fig. 
9). The detection fidelity clearly depends on the number den- 
sity of lenses, for which here we assume the fiducial value of 
3 arcmin"^. An increase in the lens population size will increase 
the number of galaxy pairs in the estimators for the angular 
correlation functions and, therefore, reduce measurement noise. 
Since Psg(k,x) and Pg(k,x) are most easily reconstructed, it may 
be sensible in a survey with poorly confined Pg{k,x) to utilise 
only those two to constrain Psg/Pg ^ rjb as a biasing function 
that can be pinned down with the highest confidence. 

The galaxy bias analysis here uses a flat prior with 
k(^.;t)l 5 2 and b{k,x) ^ 4. The prior has the effect that a 
highly-significant measurement of galaxy clustering /g „„, will 
put a lower limit of /g.mn/^^jgj, < fs^mn on the matter cluster- 
ing, where here i^ugh = 4. For the same reason, fso^mn imposes 
a lower limit of /g,,„„/(Z7high'"high) < fsjnn for rhigh = 2 to fs^,„„. 
Therefore, combining galaxy clustering and shear tomography 
aids the recovery of the matter power spectrum, as can be seen 
when comparing "independent" and "combined" lines in Fig. 8. 



Fig. 10. Relative change in the signal-to-noise in fsj,,,, if only shear 
tomography correlation functions with i j (redshift slices with 
no overlap) are used for the estimation (minimum variance). The plot 
is based on a fiducial survey with cr^ = 0. Different lines correspond 
to different z-bins. In total, roughly 30 percent of the signal-to-noise is 
lost. 



The improvement is mostly present for low signal-to-noise bands 
in the matter power spectrum. By varying the prior intervals of b 
and r, we found that the improvement is mainly due to the prior 
of r rather than b. 

4.6. Systematics 

The above predictions rely on the assumption that systematics 
in the measured correlation functions are negligible. In reality, 
several potential systematics - unrelated to the data reduction 
and shape measurements of galaxy images - have been identified 
and several points can be made: 

- Correlations of intrinsic (unlensed) galaxy image shapes 
("11") add to the shear-shear correlation functions ^^J/^ (e.g.. 
Hey mans et al. 2006). 

- The "GF'-effect is more intricate (Hirata & Seljak 2004), as it 
is non-locally produced by correlations between the intrinsic 
source image shape and the shear effect of its surrounding 
matter on a more distant source. 

- Another unwanted effect that could arise is the influence 
of cosmic magnification (Narayan 1989; Bartelmann & 
Schneider 2001). Magnification by matter density fluctua- 
tions in front of lenses changes the projected number den- 
sity of lenses; to the lowest order to /Cg i-> /<-g + Ak, where 

A :- 2{a - 1) and A^g(> /) oc f^" is the mean number den- 
sity of lenses with fluxes greater than /, and k is the lensing 
convergence by intervening matter in front of a lens. This 
adds an additional contribution to yJ'''(A6') of the order of 
A{k(6)k(6 + A6)} and a contribution to cJ'^H AQ) of the order 
of A^{K(9)K{e + AO)). Since A is usually of the order of unity 
(van Waerbeke 2010), both contaminations are of the order 
of the ^+(A0), and generated by matter in front of the lenses. 
In most cases, this signal is considerably weaker than the 
galaxy-galaxy lensing or lens clustering signal so that this 
eflBCt is likely to be of only minor significance. 
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Although the II- and Gl-contaminations are an obvious concern 
for this kind of analysis, the focus of this paper is the perfor- 
mance of the new estimators under ideal conditions. We there- 
fore defer a quantitative discussion of those effects to a future 
paper with the comment that the contamination can be partly 
removed by considering only slice pairs with little or no over- 
lap in the tomography analysis in ^'"I'K To achieve this, we redid 
the aforementioned analysis for a fiducial survey with no red- 
shift errors (cr^ = 0), i.e., and exactly no overlap for / j, and 
compared the relative change in signal-to-noise of an analysis 
that uses all tomography correlation functions to the signal-to- 
noise in a analysis that utilises only the / it j correlators. The 
impact of omitting that piece of information on the significance 
of fg is presented by Fig. 10: The significance drops by approxi- 
mately 30 percent at all redshifts. The contaminations can possi- 
bly be more effectively removed by means of nulling techniques 
(Joachimi & Schneider 2008). Since this is very likely to further 
degrade constraints on the 3D power spectra, we can speculate 
with the predictions above that significantly more survey area 
than 200 deg^ will then be needed for a proper reconstruction of 

4. 7. B-mode matter band power spectrum 

Finally, the estimator for the 3D matter band power spectrum, 
based on Eq. (11), presumes the absence of B-modes in the shear 
field, i.e., that shear-shear correlations vanish after a rotation of 
the source ellipticities e by (^45 ;= n/A, hence after the trans- 
formation e i-> e -'^^^e. The presence of a B-mode component 
in can potentially distort fg. Whether such a contamination 
is present can be verified in an analysis by estimating with 
all source ellipticities rotated in the data by (/145 (B-mode band 
power spectrum) and by estimating /™ with only one source 

ellipticity in the two-point correlator f^^^ rotated by 1^45 (E/B- 
mode cross-power spectrum). Hence, the inferred B-mode band 
powers fg and /™ that are inconsistent with zero indicate the 
presence and magnitude of such a contamination. We note that 
this only reveals possible B-mode issues in the data, but is by no 
means a sufficient and necessary condition for systematics in the 
data, since systematics could in principle only be present in the 
E-mode. 

To obtain an estimator for Pg{k,x) that is purely rooted in 
the E-mode components of the shear tomography correlations, it 
is also conceivable to devise the COSEBI formalism, which was 
proposed by Schneider et al. (2010). In this alternative frame- 
work, only the E-mode components of the two-point shear cor- 
relation functions are extracted. As the COSEBI values E„ are 
linear transformations of ^±{9) on a finite interval -t^i,) < 9 < d-i, 
they (a) can also be easily measured from the data and (b) must 
also be linear transformations of 



E^yfs 



(51) 



where the COSEBI statistics e!,'''* up to order n - \ . . . A^oider and 
pairs of redshift slices are assembled within the new data vector 
E. Therefore, one could use V instead of X and E instead of ^ 
in the previous estimator in Eq. (20) and the likelihood in Eq. 
(46) to implement a matter band-power spectrum reconstruction 
based on COSEBIs. Moreover, as shown by Asgari et al. (2012) 
and Eifler (201 1), to estimate cosmological parameters with cos- 
mic shear only a few COSEBIs modes are required. If this also 
turns out to be the case for the measurement of the band powers, 
this will dramatically reduce the size of the noise covariances in- 
volved, possibly simplifying their calculation or estimation from 



the data. It is unclear at this point, however, up to which order 
Noida COSEBIs are required to comprise essentially all informa- 
tion on Pg{k,x)- 
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Appendix A: Numerical evaluation of basis 
functions 

To calculate the contribution of a power band to a two-point cor- 
relator A,„„{0) on the sky, we note that all three families of basis 
functions in Eqs. (16), (28) and (37) can be written in the form 



e'A„U9)^ dx8(x) j ds s Us) P'^^ \^j^,xj .(A.l) 

^" K,kix)e 

For a numerical evaluation of integrals of this type, we suggest 
approximating the inner 5-integral with the sum 



f^,nn(9,x) := J dssJ„{s)P'-yjj-^^,x 

fj , i-/+i/kCf)fl 

^int— * p 

^ Y.P^'Hhx) J dssUs) 

= 2 p''\hx)[Jn{kj^ifi.ix)e)-Jn(kjMx)0) 

j=0 

where 

{xJi (x) , for n - Q 

2-2Jq(x)- xJi(x) ,for« = 2 (a.3) 

4 + (x - f ) /i (x) - 872(jf) , for « = 4 

and kj :- {kj + kj+i)/2, kj :- k„, + jAk and :- (A;„,+i - 
^m)/MntforMnt ~ 100 integration intervals. Where f*''^ (A:, ;^f) - 1 
is a reference power spectrum, the foregoing sum is exact for 
Nint - 1- The outer ;^-integral 



(A.2) 



9'aU0) = dxg(x)A„„(e,x) 



(A.4) 



can be computed with Romberg's method (Press et al. 1992). 



Appendix B: Generation of mock catalogues 

Here we outline an algorithm for making Monte Carlo realisa- 
tions of shear tomography data, including galaxy clustering and 
galaxy-galaxy lensing. 

The convergence fields are simulated as Gaussian random 
fields, whereas galaxy clustering obeys log-normal statistics 
(Coles & Jones 1991). Log-normal statistics, or statistics with 
a lower bound /Cg > -1 for the number density contrast Kg, is 
essential for a catalogue of mock galaxy positions. For the sim- 
ulation, we use 2048 x 2048 pixeP grids, spanning an angular 
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area of 4 x 4 deg^, out of which just 2x2 deg^ subfields are 
taken after performing our realisation. The simulation comprises 
N~ = 5 redshift slices. The Monte Carlo recipe is similar to the 
approach outlined in Coles (1988), Coles (1989), and Coles & 
Jones (1991), but generalised to the simultaneous realisation of 
a set of random fields, both Gaussian and log-normal, with de- 
fined cross-correlations. As a model of the 3D dark matter power 
spectrum, the prescription of Smith et al. (2003) is employed. 

For the generation of correlated random fields, we follow the 
approach of Simon et al. (2004), but now also include the clus- 
tering of galaxies as additional degree of freedom. A realisa- 
tion is performed on grids of angular size A, namely one lensing 
convergence grid k^'\0) and one galaxy number density /Cg'(0) 
grid for each of the i - I . . .N~ galaxy subsamples. The Fourier 
modes for identical £ vectors of all 2N~ grids are compiled as 
one vector 

m = (^'"(^), ~K^^\o, ■ . ■ , 4' V), ■ . ■)' (B. 1) 

and are randomly generated independently from aU other modes 
by (Sect. IV.B. in Hu & Keeton 2002) 



m 



1 



V2I 



Ls (g + ig') 



(B.2) 



where gi,g'j M{0, 1) are vectors of normally distributed ran- 
dom numbers, with mean zero and unity variance, having the 
same number of elements as v{{). The matrix L5 denotes the 
Cholesky decomposition, S = Ls , of the Fourier mode co- 
variance matrix 



P.(^) P.g(^) 
P.g(^) Pg(^) 



(B.3) 



This matrix S consists of three sub-matrices, describing (a) the 
auto- and cross-correlations between the convergence fields by 



[P.(0],7 = 



h(x) 



,X\ , (B.4) 



(b) the lenses number density and convergence cross- 
correlations 



[P.g(0],7 = 



3//,iQn 



2c2 



^ hixMx) 'H/kCr)' 



X\ ,(B.5) 



and (c) the auto- and cross-correlations between the lens number 
density fields 



[Pg(^)],7 = Pfii) , 



which needs further discussion below. For the scope of this pa- 
per, lenses are just random subsets of the source samples, i.e. 
Px^Oc) - I'x'ix)- Moreover, galaxies are assumed to be unbiased, 
i.e., b{k,x) = r{k,x) = 1. 

Case (c) in Eq. (B.3) requires a special treatment to ensure 

that a log-normal random field can be reaUsed. We obtain Pg''\{) 
from 



(0 



Jo 



fcix) 



hixY 



X 



^'\_qf(x)q^f(x), 



in three steps: 

step 1 : c^^'j\6) = ^ J^°° d£ £ JomPs'\£) , 
step 2 : w<'^'(6») = In (co'^'j^O) + 1) , 
step 3 : P'-l^\{) = In d00Jo{WW^\O) ■ 



(B.7) 



(B.8) 
(B.9) 
(B.IO) 



This transformation rule is found because a log-normal random 
field with lower limit K\n > -\ can be generated from a random 
Gaussian field /Cg via (Coles & Jones 1991) 



- 1 , 



(B.ll) 



if the two-point correlation function of the underlying Gaussian 
process is 



g(0) = U(O')k0 + ff)) = ln(l + uxM) 



(B.12) 



compared to the two-point correlations (DiniS) in the log-normal 
random field, where crl is the variance in k„. By the same argu- 
ment, however, the cross-power between galaxy number density 
and convergence grid in PKg{() remains unchanged because 



(B.13) 



for 

m,S2)j = j d6id62P{6u62)f(6uS2) , 



)\a\ + 5\(t\ 



(B.14) 
(B.15) 



exp 



[ 5\(T\+5\(T\-2u)(e)5x82'^ 



2cr](ry\-m^) 



5\Y crl 



6V1 and i/r(6») 



(jj(ff) I {0-10-2), and 



where o-\ - ww, ^2 ~ \"2l 
Pi^Xj^i) is a bivariate Gaussian p.d.f. 

After all ^-modes v({) have been realised, the grids are 
Fourier transformed to real space to yield the convergence fields 
and the galaxy number density contrast fields on a grid. The 
galaxy number density contrasts represent only the underlying 
Gaussian fields, ATg, so far and hence still need to be transformed 
into the log-normal density contrasts. This is done by applying 
the additional mapping /cin = g%-'^:/2 _ i fgj. gyery pixel. We 
acquire o"^ from the given realisation of the Gaussian random 
field. 

The mapping, however, is only an approximation because /Cg 
is actually a realisation of a smoothed Gaussian field, which is 
denoted here by an overline: 



(B-6) ^„(6I) = J d^& K„^(e')W{e-0') 



(B.16) 



The function W(0) represents the pixel smoothing function, 
which is assumed to be normalised to unity by definition, 
i.e. ^ d^6W{9) - 1. Applying the exponential mapping to 
thus assumes that a similar relation between smoothed underly- 
ing Gaussian and smoothed log-normal field ^in approximately 
holds, namely 



^ln(0) 



*„(6l)-cr;/2 



- 1 . 



(B.17) 



The relation between the smoothed fields would be exact, if the 
mapping were linear The approximation has higher accuracy 
for less fluctuating ^g and vice versa. For this paper, we em- 
ploy mock data in which galaxies perfectly trace the dark matter 
distribution at all redshifts. As dark matter fluctuations decline 
towards higher redshift, the above approximation provides a gain 
in accuracy for higher redshift bins. With the outlined method. 
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Fig. B.l. Relative bias APm/Pjif (percent) in the power spectrum of the 
log-normal random field realisation for different redshift bins as func- 
tion of angular mode wave number (. 



we achieve an accuracy in the log-normal field power of better 
than 10% and on /'-average several percent for the fiducial sur- 
vey employed (Fig. B.l). We consider this to be suflicient for our 
purposes, although a better treatment of this bias may be desir- 
able for future applications. 

Finally, the grids k^^ and serve as a basis for the ith red- 
shift slice of the mock survey: a candidate galaxy position is 
drawn randomly as the position on the grid from a uniform PDF. 
The position is accepted, if jc < 1 -H /fl''(0) where x e [0, 1] is also 
a random number from a uniform PDF, in accordance with Coles 
& Jones (1991) for grids large enough to accommodate mostly 
no more than one galaxy per pixel. This algorithm thus utilises 
a Poisson process to sample the galaxy number density contrast 
Kg^ with a discrete set of galaxy positions. As a complex source 
ellipticity, we adopt y^'^O) + e^, where y''' is the shear field cor- 
responding to (Kaiser & Squires 1993) and the random in- 
trinsic shape of the galaxy drawn from a statistically independent 
PDF. This process is repeated until the desired number density 
of galaxies for the ith redshift bin is reached. 
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